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Abstract 

In the field of holomorphic dynamics in one complex variable, hedge- 
hog is the local invariant set arising about a Cremer point and endowed 
with a very complicate shape as well as relating to very weak numerical 
conditions. We give a solution to the open problem of its digital vi- 
sualization, featuring either a time saving approach and a far-reaching 
insight. 

1 Introduction 

Let C be the complex plane. Since oo cannot be handled like an ordinary 
point here, as deserved by many complex functions, this lack is fulfilled by 
the compactification Coo = C U {oo}, namely the extended complex plane; 
anyway the representation of the neighborhood of oo is still impracticable 
here. Riemann cracked the problem by a stereographic projection of Coo 
onto a sphere of radius 1: this is the Riemann sphere C. 

Holomorphic dynamics collect the studies on the iterates f n of the func- 
tion / : C^, —¥ of given type (entire, meromorphic, transcendental, . . . ) 
and in one or several complex variables (depending on v > 1). The map / 
is of finite degree. 

We will deal with the case of one variable, where v = 1. The questions 
in this field may rise to high degrees of complication and many ones deserve 
a multilateral attack rooting into Complex Analysis, Topology, Theory of 
Numbers, Uniformization Theory. One of the most important goals is the 
study of elements not changing under iteration: the so-called invariants, 
showing up to dimension 2 at most: they might be points (dim 0), lines (dim 
1) or surfaces (dim 2). The collection of invariants with a same property is 
the invariant set. 

It is convenient to split results into the l locaV and into the 'global' branch, 
in order to have a good picture of the whole corpus. A local study investi- 
gates on the properties which hold up to a finite distance from the invariant 
set, thus inside a bounded domain; while a global approach wants those 
properties enjoyed by points being even at infinite distance, thus all over 
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the C. Branches are not disjoint and related concepts act in mutual coop- 
eration: in fact, either locally or globally, the fate of (inverse and forward 1 ) 
orbits closely relates to the nature of the fixed point 5 and its neighborhood. 

Results rely on the study of the orbits, i.e. the set of points generated 
by iterating the given function /, 

fl = /> h = /[/]> • • • > fn = f[fn-l]- 

approaching to limit sets which may consist of finitely (dim 0) , or of infinitely 
and uncountably many (dim 1 and 2) points. For such sets finitely many 
points of order k > 1, one speaks of limit cycle of periodic points S n : 

f(Si) = 52, f(6 2 ) = 5 3 ,..., /(^-i) = ft. (1.1) 

If the period is k > 1, then 5±<i<k = fk{5\<i<k)- If k = 1, we have the limit 
fixed point 5 and the expression Ql.l|) boils down to: 

6 = f(5). (1.2) 

so the fixed point can be re- framed as a cycle of period 1. 

Other cycles may belong to invariants J of dimension (infinitely many 
points), or 1 and, exceptionally, of dimension 2 2 (both, uncountably many). 

In the economy of dynamics over C, cycles of finitely or of infinitely 
and of uncountably many points play different roles which are explained 
locally in the former case, while the latter cycles are object of the global 
investigation. 

In the second case, we speak of 1 Julia sets' J: when they include infinitely 
many points, their topology is totally disconnected; when uncountably many, 
they are continuous (Jordan curves or not). 

With the caution to the casuistry of related local dynamics, a same limit 
cycle of finitely many points groups all orbits converging to it into macro- 
sets, said Hhe basins B of attraction' , also defined l Fatou sets', in honor of 
Pierre Fatou (1878-1929), who co-pioneered these investigations in the same 
times (1918-20) and independently from Gaston Julia (1893-1978), who is 
credited as the first official discoverer 3 of sets J in 1918. While the role 
of the basins is to include all those orbits to a same limit cycle, J is the 
boundary between the basins. 

Julia sets are well-known objects at all levels today, sparkling the imag- 
ination a wide range of people, from artists to mathematicians. The details 
of their suggestive and very complicated fractal shapes were disclosed to 
human eyes by the early computer experiments during late 1970s: machines 
revealed to be the indispensable aid for overcoming the long standing barri- 
ers of hand-written, rough plots available to those ancient mathematicians. 
The technologic run is continuously opening to finer results, thanks to higher 

1 Whether the iterative index n is negative or positive respectively. 
2 In this latter case, they span all over the Riemann sphere: J = Coo- 
3 (Related historical events are not fair.) Each basin relates to a same limit cycle. A 
fuller historiographic investigation will appear in 0. 
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screen resolutions and faster CPU saving from long time consuming compu- 
tations, as required to render these images. Anyway this side role shows 
that Technology is not a priority and does not fully belong to the road-map 
of this field: it is just a good tool to develop those methods mastering both 
analytical and geometrical relations - which are most wanted today -, to- 
gether with more accurate numerical precision, especially for the question 
we are going to deal with. 

Even at the graphical level, we differ local from global methods, whether 
they focus on the previous limit sets related to the two branches. 

After some introductory theory culminating into the presentation of 
what the invariant sets 7i, said ' hedgehogs', are meant in holomorphic dy- 
namics, we will illustrate the related problem by showing how most available 
techniques fail to display 7i adequately; finally, we will discuss how to solve 
this question and to code an approach via pseudo-C ++ language code. 

2 Basic theory 
2.1 The four cases 

We are going to sketch out the mathematical terms of our local environment: 
here complex dynamics are mostly interested in the orbits behavior induced 
by iterates f n (z) near limit cycles of the map / and focuses on cycles of 
order 1, i.e. fixed points 5, f(5) = 5. For an easier approach, we will assume 
fixed points from now on. Iterates are operators defined as follows in the 
forward sense (n > 0): 

f (z) = z, h{z) = f(z), f 2 (z) = /[/(*)], ... , f n {z) = f[f n -i(z)}; 

or backwards by a composition of inverse maps (n < 0): 

The starting point z of an orbit is said the seed. 

The classification of 5 is achieved by computing the modulus A of the 
first derivative /' at <5, A = |/'(<5)| and it is essential to understand the 
nature of local invariant sets, grouped here into four main classes: 

1. super- attracting fixed point, when |A| = 0; 

2. attracting fixed point, when < |A| < 1; 

3. indifferent or neutral fixed point, when |A| = 1 = |e 27rj6l |; 

4. repelling fixed point, when |A| > 0; 

Reader's understanding can be lessened by splitting this basic classifica- 
tion into two extremal cases, (super-)attracting and repelling entries, their 
directions are opposite but dynamical properties are very similar. With it, 
one understand that the indifferent case stands at the middle way and it is 
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the conjunction point among those ones, either in a inclusive way by pre- 
senting subcases where all dynamical features are shown at the same time 
and in an exclusive way, i.e. showing features which are not enjoyed by both 
extremal cases. 

In fact (super-) attracting fixed points can be reached by forward iterates, 
whereas Julia sets include cycles of repelling points and reached by backward 
iterates. Obviously, the dynamical characters (contraction, repulsion) of 
entries andEJ hold for the direct function f(z), while they are inverted 
for f~ 1 (z). Indifferent points deserve a separate discussion, where even 
the concept of 'limit' shall be carefully applied, because not conventionally 
intended like in the other cases. 

While entries , HJ) and^J) just need the modulus value, case|3J) requires 
a thorough investigation: the modulus is insufficient to distinguish the more 
complex casuistry here, as illustrated in table 12.1 I another parameter is 
demanded and our next (and immediate) choice is the angle < 9 < 2tt, 9 € 
R. Its numerical properties rule out the dynamical characters of neighboring 
orbits about indifferent points. The main separation, into the two great 
classes of 9 € Q (rationally indifferent 5, parabolic case) and of 9 € M\Q 
(irrationally indifferent 5, elliptic case), is followed by a number of sub- 
cases. The former relates to one only invariant set, namely the Fatou-Leau 
flower (see def. at p. 12%]) . and does not branch out; while 9 € R\Q generates 
a richer variety whose local dynamics get far more complicate: here a second 
level opens to local invariant sets, distinguishing for the numerical properties 
of 9, some of which are extremely weak and shunning the machine finite 
digits computation. The goal of this work is understand how and if such 
latter cases might be actually attackable in particular graphical terms which 
evince their local dynamics. 

2.2 Analogy with linear models 

One crucial tool to study the local dynamics is the Schroder functional 
equation (SFE): 

^[f{z)]=aty{z)\, (2.1) 

where ip(z) is an invertible map. Without loss of generalization, let the 
origin be fixed for /. If we replace a(z) with A = f'{z) - the meaning of the 
first derivative strengthens the application of SFE to local problems, then 
(|2.1|) turns into this version: 



yj[f(z)] = AV(z). (2.2) 

Given a sufficiently close neighborhood D of S, the left diagram below, 
as well as the right one extending to iterates, commutes when SFE holds: 



D -> 


f(D) 






C ^ 


c 



D -> 


r(D) 






c — -> 


V 
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The major goal of (|2.'2[) is to set up an analogy between the behavior of 
iterates f n (z) about 5 and an easier model; so, if SFE commutes, there is 
a local change of coordinates where /(0) = 0,f'(z) = A: the local behavior 
of orbits f n (z) can be studied by means of the simpler investigation of the 
linear map X n z. 

INDIFFERENT FIXED POINT/CYCLE 

A = crm, iai = i 



A = e 



RATIONAL 
27rifl n r- 



FATOU-LEAU FLOWER 



CENTER-PROBLEM 



IRRATIONAL 



A = e 2 



6R\( 



LINEARIZABLE 



SEMI-LINEAR IZ ABLE 



NON-LINE ARIZ ABLE 



SIEGEL DISC 



HERMAN RING 



Simply connected case 
Doubly connected case 



X 



HEDGEHOGS 



Table 2.1 : Resuming the neutrality. A diagram illustrating the classification 
of all dynamics for indifferent fixed points. 

Local theory had proven that SFE holds only for the (super-) attracting 
and the repelling case, together with the rationally indifferent and part 
of the irrationally indifferent cases (see list at p. GJ). If the local change 
of coordinates into Xz applies, D is a Bottcher or a Koenigs' domain, 
depending on A = and A < 1 respectively. For A > 1, local dynamics are 
repelling and, locally speaking exclusively, this case can be regarded as the 
converse of the attracting one; anyway here Fatou and Julia showed that it 
makes more sense 4 of studying repelling orbits and the related invariant set 
(as resumed in diagram 12.2 (1 in global terms: this is the Julia set, defined 
as the closure of all repelling cycles and being the common frontier of all 
basins B of convergence. 

The discussion of |A| = 1 requires a much more delicate analysis: first, 
its re-writing in the formal terms of ()2.2|) is 

iP[f(z)} = e 2me iP(z) => ^[f n {z)\ = e 2 ™ en iP(z), (2.3) 

where the version on the right refers to iterates f n (z), conjugated locally to 
the family of complex rotations e 27T8m z. Problems arose when it was seen 
that (|2.3|) does not always commute, given 9 € M, whether it does or not is 
up to the numerical nature of 9. For example, it never commutes for 9 £ Q, 
while it does when 9 £ R\Q satisfies this Diophantine condition: given two 
numbers r > 0, k > 2, 

\9-p/q\>r/q k , (2.4) 

4 According to their researches, one evinces the boundary role of repelling fixed points 
and cycles in the global dynamics all over the Riemann sphere. The former researchers 
did not score this goal, because they regarded repulsion like the converse of attraction. 
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Table 2.2 : Classification by motions. This is the classification of 
invariant sets for holomorphic dynamics in one complex variable by the 
three elementary motions in the plane. 



for every rational number p/q where p, q € N. Then f n (z) are said to be lin- 
earizable into a complex rigid rotation or, simply, ' linear izable'; topologically 
speaking, there exists a neighborhood of 5 which is conformally isomorphic 
to a disc and where the dynamics of f n {z) are rotational, the Siegel disc. 
When (|2.4j) holds, 6 is poorly approximated by rational numbers; otherwise 
it is goodly approximated. It is clear that rational numbers (parabolic sub- 
case, 6 € Q) are not Diophantine, together with a subclass of irrationals 9, 
the Liouville numbers - not satisfying (|2.4[) . In both cases, SFE fails and 
quits to be useful for investigations, whose achievement shall necessarily rely 
on other tools. Therefore one summarizes the previous concepts as follows: 

• when SFE commutes in a sufficiently close neighborhood of (super- 
) attracting, indifferent and repelling points, the local behavior of iterates 
links to only one among the three elementary motions (either contraction, 
rotation, repulsion respectively) [ 5 ] with regard to the fixed point 6; elemen- 
tary motions also include the translation f(z) : z + k and the identity map 
f(z) : z, known to apply at the fixed points. 

• so SFE classifies orbits dynamics by contraction, rotation, repulsion 
exclusively; it turns out that, as SFE fails, local dynamics do not involve 
only one elementary motions, as explained in diagram 12.2 I then Fatou- 
Leau flowers and Hedgehogs consist of compositions of more than one 
elementary motion. Hence the failure of SFE shuts the door to linearizability. 



The fourth linear motion, translation, refers back to contraction or repulsion cases, 
depending on the vector direction approaching or leaving S respectively. 
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2.3 Hedgehogs, the next stop 

Stepping in the second sub-level of table 12.1 1 we assume to deal with holo- 
morphic quadratic germ 6 

f(z) : e 2m6 z + z 2 + ... . (2.5) 

Iterates f n (z) can be thus expressed in this form: 

f n {z) : e 2m9n z + 0(z k ), n > 1, k > 2 (2.6) 

If compared to (|2.3j) , the question boils down to the study of the 9 values 
so that 0(z k ) vanishes or not, together with the related topological dynam- 
ics. It is straight-forward that 0(z k ) = allows to conjugate locally the 
iterations to the rigid and aperiodic rotation f n (z) : e 2wlSn z (fig. 12.3 ]/ A) in 
a sufficiently close neighborhood of 5. Topologically speaking, if 9 enjoys the 
diophantine condition l|2.4|) . the invariant set is a disc-shaped neighborhood, 
the Siegel disc S centered at 5, which Julia termed center like Poincare did 
for the analogous points in the field of differential equations: 

Definition 2.1 (Siegel Disc) Let f be a non-linear complex function (12.5(1 
and let 5 be an irrationally indifferent fixed point, so that \ f'(5)\ = e 2m ® . 
If 9 enjoys the diophantine condition (|2.4|) . the SFE (|2.3|) commutes and 
there exists a sufficiently close neighborhood S of 5, so that S is a simply 
connected component of the basin and analytically conjugated to an aperiodic 
rotation e 2lTie of the unit disc D. ([IT], p. 117) 

What happens if this condition is no longer met by 9 € R\Q? Recent 
Mathematics has been stressing that the concept of l Natura non facit saltus 1 " 1 
rules out for dynamical systems. Researches ( |3 El E3 EH- to quote a few) 
showed the existence of irrationals 9 dropping gradually off the diophantine 
condition (|2.4|) : the Siegel disc S slowly turns into a new topological family 
of invariant sets, namely the hedgehogs Ti. We cannot expect immediate 
changes in the local geometry: S squeezes into a smaller disc and neighboring 
orbits inside B\S get very complicated; if S is not maximal, we say that Ti. 
has a small linearization area. In fact, with regard of the angle value 9 
involved, S may shrink to the fixed point 5: here Ti has no linearization 
area and 5 is defined a Cremer point 8 . 

During 1920s and 30s, Cremer showed [7] that the set R\Q of all irra- 
tionals can be partitioned into two disjoint subsets, Ai,J\f: 

• let Ai be the set of angles 9 so that T> c is enjoyed, then a maximal 
Siegel disc S exists; M. is a full Lebesgue's measure set; 

• let C ^ N be angles whose value is a Liouville number. TV is a null- 
measure set. 

6 The current result apply to this hypothesis. For higher degree germs, they are sub- 
jected to changes. 

7 Ancient Roman phrasing by Linneus, alluding to Nature evolution by gradual and 
infinitesimal changes. 

8 In honor of Hubert Cremer (1897-1983) who first proved their existence when 8 does 
not meet the Diophantine condition. 
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Figure 2.1 : Singular dynamics. Tracking down one orbit inside the 
basin of attraction with a Siegel disc S. The red and blue points mark 
the seed z and the final iterated image point z n respectively. This figure 
was plot by nested circles equi-potentials. 



In the course of his proof on the existence of non-centers for rational maps 
[7j, Cremer also showed that the iterates R n (z) give rise to periodic cycles 
of growing order n in the neighborhood of the non-linearizable indifferent 5. 
He also found that they accumulate at 5 while 

liminf a ^/\f3 n - 1| = (2.7) 

n=l,2,... 

which holds for 

liminf \(3 n - 1| = and liminf (T = 1,0 = e 2nie . (2.8) 

n=l,2,... n=l,2,... 




(A) (B) 

Figure 2.2 : Wedging the bounded basin (I) Two hand-made drawings of 
the hedgehog entering the bounded basin (shaded region). The wedging white 
region is the basin of the point at infinity. 

Since the number of (super-) attracting cycles is finite (according to Koenigs' 
conjecture in 1883 and proved to hold by Fatou and Julia) and the Julia set 
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J nature (closure of repelling cycles), most of the periodic points in these 
cycles have to be repelling 9 so that, roughly speaking, [5] their accumulation 
shows that, as n grows, J wedges (see figs. 12.4 I and 15.8 \ the bounded basin 
of attraction: if n has no upper bound, then 5 € J, otherwise the wedging 
action of J stops at a certain distance from 5 for n — > K < +oo. The zero 
limit in (|2.7|) is reached when a is rational 10 or it is a Liouville number; 
hence a poorly rational approximation - for the diophantine condition for 8 
- implies rotation around 5, because cycles do not accumulate at it. Inspired 




(A) (B) (C) (D) 

Figure 2.3 : The ideal journey. Simple figures of the local invariant sets related 
to the numerical properties of the angle for indifferent points: from the Siegel 
disc (A) to a Fatou-Leau flower (D). 

by continuity, it is suggestive to think of arranging an ideal journey, along the 
number properties so to visit all invariant sets arising for indifferent fixed 
point, from (A) to (D), that is, numerically speaking in terms of 9, from 
M\(Q> to Q. During the journey, the Siegel disc S squeezes (fig. 12.3 ^ B) until 
it disappears (fig. 12.3 ^ C): for both cases B and C, the invariant sets are 
'hedgehogs' TL, when S is not maximal or has empty interior. This additional 
look-up at the journey is inspired by the strong unifying power offered by 
hedgehogs theory, which links apparently far dynamical configurations. The 
mathematical definition of the hedgehog TL is |14j as follows: 

Definition 2.2 (Hedgehog) Given a neighborhood U of an irrationally in- 
different point 5, so that the holomorphic map f is univalent on U , the 
hedgehog TLisaninvariantcompactum 11 , so that f is not linearizable or has 
linearization domain relatively compact in U . 

For (|2.5[) . the current research investigates on the semi-local dynamics 
about TL and thanks to the results by Perez-Marco in 1990 ([Hj, P- 123) 
one knows that TL may show up in these different three characterizations: 

1. hedgehogs TL are locally linearizable and with no small cycles; 

2. hedgehogs TL are not locally linearizable and with small cycles; 

3. hedgehogs TL are not locally linearizable and without small cycles; 

9 Cremer never questioned on the nature of the accumulating cycles at <5 throughout 
his works. 

10 This case is out of our context. Anyway the Fatou-Leau flower dynamics are induc- 
tively helpful to watch how the basin to oo wedges the other basin, so that the Julia sets 
attaches to <5, assumed the germ f(z) : e 2 ^ 16 z + z n , where n £ N + ,n — > oo and 8 G Q. 

11 A compactum (plural, compacta) is a compact metric space. 
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Figure 2.4 : Hedgehogs generation. An illustration of 3 
steps (different shades of grey) of repelling cycles accumulation. 
The starred-shape, due to the cycles getting near 8, as well as 
wedging action of the Julia set J, are evident. 



The ' small cycles' property refers to the possible existence of infinitely 
many cyclic orbits in a sufficiently small neighborhood of the irrationally in- 
different 5. It is said that l 5 has the small cycle property 1 or l is approximated 
by small cycles 1 . Other results on small cycles are included in |19U2U| . 



3 Entering the computer graphics 

3.1 Features of ordinary methods 

Most available methods are global and devoted to the display of Julia sets J. 
One first notices that they do not often 12 focus on the analytical properties 
of J. Both its location and shape come out from the human eye perception 
of the different colors distribution inside a close neighborhood of J. In 
mathematical terms, these methods do not draw J: they paint the Fatou 
set T = C\J. Therefore they are not detective, but deductive: focusing on 
the complement set they deduce J from it. 

3.2 Palette of colors 

Experience induces this consideration for color figures: the (ascendent or 
descendent) ordered palette of shades 13 fits for understanding the conver- 
gent or divergent 14 movement of iterates between distant neighborhoods of 
the fixed point, because the shades sequences play as graphical indicator 
of the dynamics, but they fail to evince the local dynamics: the conver- 

12 The approach by inverse maps achieves it, but it is not exportable to any maps because 
inverse ones cannot be always retrieved. 

13 For example, sorted by gradient: shades of blue or tones from green to violet. 
14 With regard to the nature of 6. 
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gence/divergence rate about the indifferent 5 is very slow and iterated do- 
mains match to very close colors in the palette, looking like the same. 

3.3 Convergence criteria 

The fate of iterated orbits is tested at each f n (z) : z n , for preventing infinite 
looping. This is commonly achieved by a trapping disc T> of given radius r; 
then \z n \ < r is tested. This approach, aiming to understand when (in terms 
of the value of the iterative index i) orbits might escape T>, is commonly 
known as 1 escape time'' algorithm. Some pseudo-code follows: 

#include " complex. h" // this is a class handling complex numbers 
// downloadable from author's site 

complex z, next_z ; 

complex c(0, -1); // this is the complex parameter z = 0.0 - l.Oi 

t = 50 ; // top iteration index to prevent infinite looping 

r = 2.0 ; // the radius of the trapping disc 

for( int i = 0; i < t ; i++ ) 
{ 

next_z =z*z+c; //we assumed the quadratic iterator 
if ( abs( next_z ) > r ) break; 

z = next_z ; 

} 

take-some-color-value-f rom-the-iterative-index-or-f rom-the-point-z ; 
plot-z-on-the-screen ; 

Another approach relies on the convergence properties in the neighbor- 
hood of 5 and uses the distance between two successive iterates as test 
condition: 

\fi(z) - f i+1 (z)\ <e, ,e>0 (3.1) 
It is a variant of escape time and defined the 1, approximation 1 method: 

#include " complex. h" // this is a class handling complex numbers 
// downloadable from author's site 

complex z, next_z ; 

complex c(0, -1); // this is the complex parameter z = 0.0 - l.Oi 

t = 50 ; // top iteration index to prevent infinite looping 

e = 0.00001 ; // the distance ranging between < e < 1 

for( int i = 0; i < t ; i++ ) 
■C 

next_z = (2*z*z*z+l) /(3*z*z) ; // this is the transformed map 

// of z~3-l by Newton's method 

if ( abs( next_z - z ) < e ) break; 
z = next_z ; 

} 
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take-some-color-value-f rom-the-iterative-index-or-f rom-the-point-z ; 
plot-z-on-the-screen ; 



3.4 Obsolescence with local invariant sets 



We are going to discuss ways to bypass or lessen the two problems in subsec- 
tions 13.21 and 13. 31 but even any improvement in such two directions would 
be still insufficient for our goals. First we deal with colors: sequences shades 
showed up to unfit, so might a randomly generated palette help? In the 
color figures of table 13.51 we iterated the neighborhoods about fixed points 
of different nature and painted them by the random palette. For figs. 13.5 | /A 
and B, we picked up the square map f(z) : z 2 , having a super-attracting 
fixed point 5 at 0; in (C), the Newton- Raphson method was applied to 



1, having 3 attracting points on the unit circle <9D. One no- 




(A) 



(B) 



(C) 



Figure 3.5 : Attempts of improvements. The random palette coloring for 
attracting fixed points. More than the shading gradient approach, the attracting 
dynamics are evident here from the nested discs shrinking up to points. 



tices that this works finely with the neighboring dynamics, whose shrinking 
behavior to the attracting fixed point(s) is displayed step by step, or analo- 
gously in the neighborhood of a repelling point. 

In the figures of tables 13.6 I and 13.7 1 we painted all invariant sets (refer 
to fig. 12.3 I which could possibly arise in the neighborhood of indifferent 
points 5 by ordinary method, with emphasis to colors sequence and values 
accuracy. 

The iterations of f(z) : z + z 4 yield a 'Fatou-Leau flower' with 3 petals 
meeting all together at the origin. Either with the aid of random palettes 
or as infinitesimal values of e are set into inequality (|3.1|) . one cannot evince 
the local dynamics, but only the basins shape. And no benefit is drawn 
from the experiments with irrationally indifferent points. In fig. 13.7 | /A and 
B, the random palette and a very high number of iterates were respectively 
tested to work with a Siegel disc case (it should appear in the green basin). 
The plot of an hedgehog was tested in (C, D). In particular, (B) and (D) use 
another coloring method which sets a one-to-one map between each point 
and the RGB cube, so that any complex point is univocally colored with 
regard to his location. Even if the iterative index i increases to very huge 
values, orbits cannot reach too close to 6 so we cannot get better pictures 
of it. 
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(A) (B) (C) (D) 

Figure 3.6 : Ordinary attacks to rationally ... . Different methods for the 
Fatou-Leau flower. 




(A) (B) (C) (D) 

Figure 3.7 : . . . and irrationally neutral points. Dynamics near a Siegel disc 
(A), (B) and an hedgehog (C), (D). 



4 Off to 'composite dynamics' 
4.1 The core questions 

What does to draw a local invariant set mean? And, in a larger sense, what 
do we want to get while drawing it? 

Without basing upon good responses to these questions, it does not make 
so much sense to go on with our attempt. As remarked in table f2T2~l at p. 
H3 holomorphic dynamics have different invariant sets. 

Julia sets J are easy to display: the only issues are merely technical 
(methods, palette colors, ...). In local terms, the neighboring dynamics 
about points of J can be easily collected under one easy concept: orbits, 
whose character is initially repelling near J, run across the basins B of 
attraction up to a (super-) attracting cycle (here, a fixed point 8). 

Some local invariant sets I s around 5 enjoy more or less complicated 
dynamics. The table 12.2 I offered an overview of cases, where we noticed 
invariant sets endowed with 'composite dynamics', that is, involving more 
than one elementary motion: this is a bunch of information which is harder 
to be pulled out than when SFE commutes and there is only one local 
elementary motion to deal with. One should also take in account that the 
convergence speed of neighboring orbits may get slower as they reach closer 
to 5: in this case, understanding such motions via colors is impracticable. 
Here one needs to evince either the (1) boundary dZ s (retrieving the 'shape') 
and the (2) inner motions (figuring out the 'behavior' of orbits). 

Reaching (2) needs to consider that the elementary psychology of Space 
relates the perception of motion 'to lines in the principal directions of curva- 
ture, which may communicate surface shape better than lines in other direc- 
tions 1 (|H|, p. 43). This is a common feature of the field of Non-Photorealistic 
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Rendering (NPR) and associated to the concept of 'perceptually efficient im- 
ages', which is a visual representation emphasizing important features and 
minimizing superfluous details: our rendering are somehow skinny. For our 
context, we find useful to filter this concept by our needs: we want an 
analogously-shaped geometrical model of the given vector field, which is 
capable of evincing (1) and (2). 

4.2 Obstructions (I) : no ticket for the ideal journey 

The problem of drawing hedgehogs on a computer was first posed to the 
author by professor Perez-Marco during an informal meeting in summer 
2002. In general, the degeneration of a value € M\Q into 8 £ Q by finite 
digits machine computation is quite obvious. 

Assuming the formula (|2.6|) . one major issue is to check whether the 
resulting Q = 6 n may approximate irrationals poorly or not. One would like 
to distinguish the formal value 9 with the input one a in practice, where 
a ~ 9 necessarily holds: here \a — 8\ = e > 0. Analogously, one has 
\a n — 6 n \ = p > e > for iterates, thus the error magnitude grows with 
the iterative index n. Resuming, we are not contented by the value 9: how 
could we reduce p, if possible? In fact subtle issues may come up about the 
possibility: first, Cremer's partition of irrational numbers shows that values 
of the null-measure set Af are almost impossible to pick up by a machine 
with arbitrary precision, in order to guarantee the chosen invariant set: for 
example, hedgehogs 7i with no linearization area relate to extremely weak 
numerical conditions under the decimals cut-off, so that one could never 
watch them by relying on numerical computations. So would it make sense 
to approximate a number a n G R\Q via another irrational 6 n which possibly 
enjoys different numerical properties? We reply that the value a n may give 
rise to a local invariant set differing, more or less evidently, from what we 
really wanted from 9 n . We observe that from Lebesgue null-measurability, it 
does not follows that M should be pointwise and totally disconnected; rather 
it consists of disjoint intervals whose amplitude is positive but infinitesimal, 
say, of width e. If machine approximation may actually work under such 
interval e, we might think of being able to get very good approximation of 
the original value a n and of finding the related invariant set with adequate 
precision. 

Without this chance, the approximation problem is weakly attainable on 
a numerical basis; is there an alternative attack relying on the approximation 
of properties? We recall that approximation in itself refers to another entity 
[3, inexact and imitating (not copying) the original a; so, dealing with prop- 
erties, we questioned on whether it is plausible to seek one property relating 
to a third value [3 and so that 

\a — f3\ < \a — 9\ 

holds? We would enter an ad infinitum process generating an infinitely 
many values £, 7, . . . , r, and so that 

\a — (3\ < \a — £| < \a — 7] < • • • < \a — t\ < \a — 9\ 
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and therefore infinitely many properties V standing between the formal value 
a n and the practical value n and even nested as follows: 

Vp c P € c Vj c ... c V T c V e ■ 



Question 4.1 Is the number of numerical properties finite? If so, how many 
they are? 

We conjecture that the set of properties V, involved in the local lineariza- 
tion around an irrationally indifferent 5, is finite: because the topological 
configurations of hedgehogs are also finite. 

4.3 Obstructions (II): practice is what matters 

Now with regard to the practical computations involved in the iteration 
process, we can find three classes of obstructions: 

1. Statistical: due to the A4,Af sets distribution over the real interval 
[0, 1] , it much easier to pick up a 9 € M. than one Q G M. And due to 
above approximation, this task gets harder than ever. 

2. Numerical: numbers 6 £ A4 are more robust to such approximation 
attack than £ € J\f; in fact, under iterates, Liouville numbers £ tend 
to be turned into new values which are Diophantine irrationals; new 
resulting values cannot goodly approximate rational numbers and then 
one falls back into the Siegel disc case. 

3. Procedural: Liouville's condition is fundamental to plot hedgehogs, 
but it is also weak to be preserved under iteration. The movement 
rate or orbits in the hedgehog Ti gets as slower as closer to the fixed 
5. Thus one would be pushed to increase the iteration index n to get 
finer results but we saw the iteration process supports obstructions ^ 
and El 

Approximation and computer usage are two 'obstructions' we cannot 
overcome. Not having a valid solution at hand nor being successful to find 
it, the problem was left open until we occasionally came to it during March 
2006, when we planned to find a strategy for lessening the break of irra- 
tionality. It was clear that the numerical path was impracticable. 

Our problem with hedgehogs refers back to the impossibility of drawing 
local invariant sets, where slower and slower convergence rates, as iterates are 
closer to 5, escape the possibilities offered by the standard methods: in the 
next section, an analysis on the ordinary graphical methods for holomorphic 
dynamics will be given in order to motivate their exclusion from the run for 
displaying local invariant sets and, most important, hedgehogs accurately. 



16 



Digital visualization of hedgehogs 



4.4 Quality vs. quantity 

The methods described in section |3] are structurally weak for our purposes: 
the reader shall know that they have been developed (read, customized) for 
iterates of polynomial maps, in primis the quadratic f(z) : z 2 + c, where 
z,c£C and c is a parameter. As we showed here and even for quaternionic 
Julia sets, as discussed by the author in |16U17) . they might not work finely 
or even they may fail completely if exported to dynamical systems which 
are distinct from their original context. We remark that, mostly due to the 
very slow convergence rate of orbits near the indifferent S, the escape time 
approach is inconvenient to display the hedgehogs, either if different colors 
sets are applied or if the index i is hugely increased. Ordinary methods work 
quantitatively and they fail with dynamics wanting the sharpest numerical 
precision; in addition, the geometry of these invariant sets (both in local 
and global terms) is exclusively retrieved by the value of the last iterated 
point z m , where m is any largest iterative index i available to the machine 
architecture. We may set m to the largest value then, but computations 
would be excessively time consuming. 

In the spirit of automatism helping to save time, we state that a radi- 
cally different method is required: thus a way out is not to reach to those 
neighborhoods, but to 1 be-already-there-and-work' . Even in this new direc- 
tion, the method customization is required: in fact, we judged useful to 
adopt an half-way strategy which, besides the inevitable finite digits com- 
putation, teams up with a qualitative attack by the imitation of a model 
with pre-defined shape. 

The decision on the model matured after the conclusions in section PZ"21 
where the failure of SFE for Fatou-Leau flowers and hedgehogs implied that 
such dynamics cannot be adequately described by models based upon reg- 
ular curves, such as concentric circles or sheaves of straight-lines. So we 
wanted to adopt a non-regular and, mostly, topologically equivalent model. 
Namely, this is the (holed) n* ft -branched star, as depicted in table Iq"TT~1 It is 
extremely important to remark that the following results are based upon the 
iteration of quadratic germ (|2.5[) . Actually there exist obstructions of the- 
oretical nature, stopping to extend the results from quadratic holomorphic 
germs to higher degrees. 

5 Re-elaborating the equi-potentials 
5.1 The holed n-th. branched star model 

Our approach roots to nothing new in digital graphics for dynamical systems 
and often used in complex analysis. Known under the term of ' equi-potential 
curves' 1 , the reader can find examples of it in |12j by Needham and, for 
holomorphic dynamics, in by Milnor - although it already appeared, 
but weakly, in older publications JHIE]- The further pseudo-code, reported 
to support the concept rather than being implemented as it is, consists of 
a main routine assuming the whole input finite subregion of C as a grid of 
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seeds z and filtering the dynamics of each iterated point z n via a sub-routine 
checking if an orbit can be possibly plot, according to the holed star model. 

#include " complex. h" // this is a class handling complex numbers 
// downloadable from author's site 

complex fxd_pt ; //we assume that we already 

fxd_pt.real =0.0; // know the location of the fixed point 

f xd_pt . imag = 0.0 ; // around which we want to plot the hedgehog 

double val = 0, tmp_val = ; // values which will be later 

// stored and compared 

BOOL blnitFlag = FALSE ; // this boolean flag will be 

//to check the two above [double] values 

// we assume that the four coordinates of the screen port where 
// the image will appear are stored in these variables; 

int top = 0, bottom = 320, left = 0, right = 200 ; 
complex z ; 

for ( x = left + 1; x < right ; x++ ) 
{ 

for ( y = top + 1; y < bottom ; y++ ) 



{ 



// suggestion: do a converse estimation map to come from screen 
// coordinates to the complex point z further by 'RescaleToMap' 
// so that you're assured that decimals approximation 
// does not make more than one screen point to the same 
// complex point z 

rescale-the-pair ( x, y ) to-the-pair ( z.real, z . imag ); 
//////////////////////////////////////////////////////////// 



// at this point you shall just iterate the function 

// for the required number of steps but without 

// test conditions. See previous code. 

// Just stop it by a iterative top index limiter. 



// We did not show any related code for the function 
// to be iterated because one needs more complicate 
// explanations. 

// Make sure you are working with a suitable function 

// retrieving an hedgehog. We suggest you to use the one 

// an holomorphic quadratic and indifferent germ where 

// the angle theta is represented by the continuous fraction 

// reported further. 



// We can say we coded a complex parser 

// to input any function even in the form; sorry 

// but the latex syntax is the best we can use to 

/ / arrange one example : 

// f(z) : e~(2\pi i\theta)z+z"2 

complex output_z = iterate-the-f unction f(z) ; 

// NOTE : if you input the identity map, you'll see 
// the holed branched star. 
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BOOL bDraw = Holed_Star( output_z, fxd_pt, doublefe val, 

doublefe tmp_val, B00L& blnitFlag ) ; 

if ( blnitFlag ) 
{ 

if ( val ! = tmp_val ) 
{ 

if ( bDraw ) pDC->SetPixel( x, y, ); 
val = tmp_val ; 

} 

} 

else blnitFlag = TRUE ; 

} 

> 

In the previous code, one sees a raster 15 scan of the screen by rows (y) 
and columns (x). As integer coordinates are taken per each screen point, 
they have been immediately turned into real and imaginary values respec- 
tively, so to have finally the complex point z. 

And now we give the code details for the function HolecLStar which 
draws the hedgehogs/holed star by equi-potentials point after point. In this 
code, we just presented monochromatic hedgehog (refer to table !5.lT~l A/B), 
the colored version requires one more function which we chose to not include 
here for code simplicity. 

BOOL Holed_Star( complex z, complex fxd_pt, double& val, 
doubleft tmp_val, B00L& blnitFlag ) 

{ 

#include " complex. h" // this is a class handling complex numbers 
// downloadable from author's site 

#define PI 3.14159265358979323846 // we need PI to be as 

// precise as possible 

complex tmp_z = z ; //we need another variable for 

// storing a temporary value later 

double branches_number = 12 ; // this explains itself ! 

double offset = 0.5 ; // this is the radius of the 

// disc hole in the star 
//if set to zero, the disc is empty 

double existence_interval = 2 . * PI ; // the values range in which 

// compute all branches location 
// in radial terms 

double potential_rate = existence_interval / branches_number ; 



15 This is a technical term in computer science to indicate that an image is regarded as 
a mesh of points distributed in rows and columns; each point is associated to a triplet of 
values. The first two values are the unique pair of coordinates which define the location in 
the mesh. The third value refers to the point color, usually defined in the RGB additive 
model. 
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tmp_z -= fxd_pt ; // performs a translation of all points 

// so that the fixed point maps to the 
// origin and the further computations 
// get easier 

double v = tmp_z . angle () ; 

double out_level = ( f . euclidean_dist ( z ) <= offset ) ? -1 : v / potential_rate ; 

int ol = (int)out_level ; //we need an integer value for the 

// output level to be eventually used 
// for color the point 

if ( IblnitFlag ) val = ol ; //if the flag was not set, 

// then this means that it is the first ever 
// value to be stored. 

else tmp_val = ol ; // Otherwise, it does not. See before 

// in the main function 

if ( ol < || ol > branches_number ) // check for values not 

return FALSE ; // ranging out of the interval 

return TRUE : 



We have now a bunch of code to scan a region 1Z C C, so that HolecLStar 
classifies the n-fold image region f n (Jfy via the star model. 



I 



1 

X 

4* 



(A) 




(B) 



(C) 



Figure 5.8 : Wedging the bounded basin (II). Three close-ups of the Julia 
set J for the function 12. 511 with 9 expressed by 15.11 , rendered by the holed star 
model. The iteration index was set at (A) 150, (B) 300 and (C) 1000 respectively. 
White disks are preimages of the hole in the star, while colours are changing due 
to iteration values attained by the forward images at different ranks. The wedging 
action of the hedgehogs generation is evident. 



5.2 Unveiling the code 

Here we explain the elaborations worked by part of the code in section 15.11 
First one shall focus on this line: 

double out_level = ( f . euclidean_dist ( z ) <= offset ) ? -1 : v/potential_rate ; 

if the point z, resulting from the iteration process, is farther from the fixed 
point f than a preset distance offset, the ratio v/potentialj:ate shall 
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be computed to know whether one paints that iterated point black or not; 
otherwise, the negative value —1 will be used as a flag, indicating that the 
point must not be tested for further painting. The variable offset stores 
the radius value for the hole in the star model. 

The value v is the angle of the last image z in the current orbit and 
helping to have that point location in relation to the branches. One also 
notices that 

double potential_rate = existence_interval / branches_number ; 

which, applied to the first line of code, let v / potentialjrate turn into 

v • branchesjiumber v 

index = = • branches_number 

exist ence.interval exist ence.interval 

the ratio factor in the right member is normalized in the unit interval [0,1]; 
when multiplied by branches_number, one finds that: 

< index < branches_number , 

so it ranges inside the interval of arbitrary input number of branches. If 
decimal values are retrieved, they shall be rounded off to the integer part, 
because decimals do not make sense in terms of branches index. For example, 
given 31 branches for an holed star, like in table 15.10 1 we expect that 
the index ranges in [0,30]. The product with the branches_number yields 
the branch index for the last computation of the n-iterated image z. It 
is reasonable that watching hedgehogs pictures one first wonders: given a 
white background, how come are black lines painted exclusively? And, in 
addition, since the black pixels refer to those only points to be painted, 
why are others not painted? One side implication from a general method 
based upon equi-potentials is that screen points are plot only if the current 
equipotential is trespassed. 

In our computational model, the geometry locus of the curves is the 
qualitative part. This locus is often detected by a formula: for example 
\z — C| — R is known to detect a circle centered at £ and if '=' replaces with 
'<', we deal with the open disc. 

In order to meet machine requirements operativity, we shall deal with 
quality in terms of quantity: therefore, each curve matches one only index 
value. When we point out to a curve or to a branch in the star model, we 
want its index essentially. With 'an equipotential level was trespassed 1 , we 
mean that the current index value has changed. For a set of points retrieving 
a same index, only the first one shall be plot; again, if the indexes of two 
consecutive points are different, the second point is plot. Algorithmically 
speaking, one stores the last index and, if it equals the previously stored 
value, no equipotential is was trespassed; otherwise, it was and one paints 
the related point black. So the implication is straight-forward: 

Different indexes — > Different equipotential — > Trespass — > Paint the point. 

For example, if we associate each screen point to its resulting index, we 
find a chain of indexes indicating how pixels and the entire screen row is 
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regarded in this raster process from left to right: 

1-1-1-2-2-2-2-2-2-3-3-3-3-4-4-4-4-2-2-2-5-5-5-5-3-3-3-3-3 

The index in bold makes the difference and arrests the previous sub-sequence 
of same indexes: the related point will be painted black. Technically, this can 
be easily understood by handling the boolean flag blnitFlag of the main 
routine at section lo.ll The reader would take care that as the blnitFlag 
is set, the values for comparison were stored in the two variables val and 
tmp_val, as shown here in the following resume of the previous code: 

BOOL bDraw = Holed_Star( output_z, f xd_pt , doubleft val, doubleft 
tmp_val, B00L& blnitFlag ) ; 

if ( blnitFlag ) 
{ 

if ( val ! = tmp_val ) 
{ 

if ( bDraw ) pDC->SetPixel( x, y , ); 
val = tmp_val ; 

} 

} 

else blnitFlag = TRUE ; 

It is clear that the point is painted black (RGB value is 0) only when 

if ( val ! = tmp_val ) 

is satisfied. If so, the old equipotential was trespassed. We jumped into 
another equipotential level (iterations are a discrete model of a continuous 
map and jumps are solely achieved onto the successive or onto the preceding 
equipotential level). And we take the new value as the referring one by val 
= tmp_val ; 

5.3 The pro(s) and con(s) 

Benefits: the fate of orbits f n (z) : z n , which finally land on a same equipo- 
tential level 1Z - the bounded region radially distributed around 5, is grouped 
by the seed points z under the same color, according to our approach to 'be- 
already-there-and-work' (p. I16|) . From the early works by Cherry [£] and the 
latest production by Perez-Marco [Tl|, hedgehogs dynamics can be simply 
resumed as follows (0, p. 33): 

"'. . . as it progressively deformed through starfish shapes [. . . ] con- 
sisting of an infinity or 'rays' emanating from O; each such 'ray' is 
a connected closed set, and 'most ' of them are arbitrarily short. A zq 
on any of these rays gives a chain (z n ) each of whose points lies on 
another of them ..." 

Hence the hedgehog boundary H of the regions 1Z wedges the bounded 
basin, which gets thinner and thinner as the iterative index n — » oo. The 
final shades distribution helps to evince the topology of Ti (figs. 15.8 \ . but a 
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Figure 5.9 : Wedging the bounded basin (III) By zooming in the red square 
region in succession, we follow the wedging action of the hedgehog inside the 
basin boundary. (A) without requiring a huge number of iterates (just 100 here) 
the contour is able to show how the hedgehog makes its way into the basin. (B- 
D) There are several smaller and smaller fjords emanating from the longest arm 
which runs towards the linearization domain (whether it is empty or not). 



major benefit comes from the contour plot of all such TZ ; contours allow 
to watch the complete hedgehog topology, by a comparatively little iterative 
index, while other methods, even if pushed to the top machine performance, 
would have not reached so far. 

Good pictures on the way TC wedges the bounded basin are granted 
because of the holed star fits the hedgehogs topology, without needing to 
input very huge indexes for having almost the same results which previous 
standard methods: in fact they suffer the very slow accumulation rate of 
neighboring orbits z n about 5. This achievement completely fills the lack of 
the approaches discussed before. 

Lacks: more than any other equipotential model, the holed star does 
not grant an 'intelligent approach', that is, capable of fitting any input 
function f(z) as best. Equi-potential curves are just as a sort of elastic sheet 
being stretched by f(z): they are used to watch the local dynamics via lines 
distribution, but an a priori knowledge about the property of given dynamics 
is required to choose the best model to depict them, among a bunch featuring 
a given equi-potential curves distribution (horizontal or vertical straight- 
lines, concentric discs, stars), each one working differently on the current 
case, thus being restrictedly indicative on the function behavior. Hence we 
may look at equi-potential models are differently styled interpreters and, like 
for any interpretative process, one must speak the same 'language' adopted, 
here of the dynamics of f(z), for best understanding: so the essence of any 
model is to play as both the proper syntax and as the fittest dictionary 
to achieve the correct 'translation'. Hedgehogs and holed star model do 
not escape this rule, so that the fine tunings can secure the most adherent 
results. This is not so difficult however, anyway the star model is less 'ready- 
to-go' than concentric circles or straight lines distributions of equi-potentials, 
because of one more parameter to tune, namely the hole radius, in order to 
properly plot the variety of hedgehogs, as listed at p. E3 



Which is the same as plotting the contour of the star itself, in equi-potential terms. 
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(C) r s = 0.10 (D) r s = 0.09 

Figure 5.10 : The radius test. Tuning the radius rs of Siegel compactum. 



5.4 Tuning the disc radius 

Being this one an empiric approach, one needs to tune some parameters to 
'grasp' the finest pictures of hedgehogs, letting their characters to match 
with the real figure as closer as possible. In figures at table 15.10 1 one 
tunes the maximal disc/hole radius rs for which the hole of our branched 
star model would look like a disc: the response from too large values (fig. 
15.10 ^ A) offers a strongly deformed curve indeed, while smaller ones may 
not evince the effective extension of the Siegel compactum. Thus one has to 
try to decrease the radius value r$ until the disc comes out. 

The author is working on a theory to retrieve sharp estimation of such 
radius metric, as well as for the Siegel disc ^Sj: this might be helpful to 
finally have very fine pictures of hedgehogs. 

5.5 Stars and equivalent classes 

As shown, orbits accumulate at their forward invariant limit sets. In the 
case of hedgehogs, whether the linearization disc is small or has zero area, 
the holed k-th. branched star is a model describing this accumulation in 
the same way as neighboring orbits do. The star model benefits are two- 
fold: while the regions imitate as best as possible such dynamics inside the 
bounded basin, the branch lines arc screening the shape of the basin 
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to oo. Since the latter and the bounded basin are complementary sets, in 
terms of converging orbits, tracking down their shapes equals to know the 
hedgehog topology. 

With regard to the star and the regions between, orbits are parted into 
equivalence classes 

£\ , £2 , • • • , £k , 

where k is the number of the branches. The basin Boo to 00 does not belong 

k 

to any class £, thus it does not belong to £ n , since the union is the 

71=1 

bounded basin B$ itself. In the holed star model, B$ is mapped to the union 
of the regions between the branches. 

The branches of the star model elongate from the outer region up to the 
disc, or to the center point if the hole has zero area: this is not merely obvious 
and plays relevantly during orbits classification. In fact, in a similar fashion 
but applied to iterates, the imitation of the branches distribution allows to 
track down the wedging action of Boo up the boundary of the hole or, again, 
to the center of the star when the hole is empty: even here, consequentially, 
the wedging action by B^ to B$ helps to deduce the hedgehog shape ( This 
is the reason why one should already know the possible existence of the hole, 
i.e. of a linearization area with positive radius.) As remarked differently in 
section 15. Ml one great benefit is the imitation up to the hole boundary. In 
fact, nearly iterates get slower and slower to converge or even their might 
not converge at all for hedgehogs. The imitation of branches distribution 
allows to break down the slow convergence barrier and the resulting line 
show the orbits fate, even after a relatively 1 slow number of iterates (see 
figs. KWl . 

In fact the equivalence classes allow to reconstruct the basins distribution 
in the same fashion of the star model: so one also understands why it also 
works finely with the Fatou-Leau flower, a similarly shaped vector field, 
where one basin wedges the other. 

5.6 Displaying the Fatou-Leau flowers 

The holed star immediately worked very finely with topologically equiva- 
lent invariant sets, so one might like to consider the benefits from this new 
graphical construction and focus on the Fatou-Leau flowers in table 15.12 I 
before moving definitely to hedgehogs. 

Definition 5.1 Fatou-Leau flower. If z is a fixed point of multiplicity 
n+1 > 2, then there exist attracting petals Vi, V2, ■ ■ ■ , V n f or the n at- 
tracting directions at z, a repelling petals V[ , V 2 , ■ ■ ■ , V n for the n repelling 
directions, so that the union of these 2n petals, together with z itself, forms 
a neighborhood Nq of z. (^2> P- 105) 



But the smaller neighborhoods to blow up are, the finer they need to be displayed 
and thus the larger iteration index shall be. 
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(A) (B) (C) (D) 

Figure 5.11 : I am the star. The (holed) 7i*' 1 -branched star in monochromatic 
and colored versions. 



One first notices that straight lines meet at the origin and they lie at the 
attracting directions, whereas the divergent ones are deformed into petals 
shaped curves. 

It is useful to remark that straight lines leading to the origin are adjacent 
to converging orbits: the imposition of the holed star model lessens the dif- 
ficulty of reaching closer neighborhoods of 5, owing to the slow convergence 
rate inside. Thus the computation and display of such flowered invariant 
set may be a good entry point for realizing the benefits of such a qualitative 
approach and the need of tuning properly the parameters of the star, such 
as the radius hole (set to for flowers) and the branches number, especially 
for hedgehogs. 



5.7 Experiments with hedgehogs 

This latter aspect should be greatly taken into account for hedgehogs too. 

(p. 123) was rather inspiring to try enhancing their graphics by 
focusing on an imitation graphical model based upon equi-potentials; this 
book also includes one formula and one figure of a quadratic holomorphic 
germ providing an hedgehog. Here the irrational angle 6 can be expressed 
in terms of the continuous fraction expansion 



' 20000+.. . 

The choice of such model was not at random: the different approaches 
in table IS"T3~1 lack of displaying the true shape, so that at last the holed star 
was chosen to best imitate the general hedgehog shape under iteration, when 
the Siegel compactum is empty or even arbitrarily small. At this concern, 
setting the right number of branches improves the figure. Our approach 
wants to exploit the equi-potentials method by looking at the deformation of 
the (holed) n tfe -branched star. One notices the existence of curves indicating 
the wedging action into wedging Bg; this was not evident from figures ETT^ C 
and D by ordinary methods. Although we cannot depict a detailed shape of 
fjords, we may at least understand where they will get to. 

As we showed, this is not an intelligent method. Thus finest drawings 
require to already know the hedgehogs features, like the the radius of the 
Siegel compactum and number of fjords, regions where penetrates Bg) 
and modelled by the star branches. 
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z + z 5 e 2nit z + z 2 , t = 3/7 

Figure 5.12 : A bunch of flowers to the star. The application of the holed 
star to different cases of Fatou-Leau's flowers. 




Figure 5.13 : Different equi-potential models. The failures to evince the 
hedgehog shape: concentric disks, vertical and horizontal lines distributions are 
displayed in figures (A), (B) and (C) respectively. Finest figures need the most 
equivalent equi-potential model. 




(D) (E) 

Table 5.3 : Tuning the number of branches. We chose for this example the function f{z) : z + z°, which yields a flower with 4 petals. The 
best performance for the holed branched star is to have some branches adjacent to both the 4 converging and the 4 attracting directions; then 
it is wise to set the branches number to multiples of 4. The other branches are useful to understand the dynamics. Without them we should 
just see straight lines, as depicted in (A). So we have 4, 8, 16, 32 and 64 branches for figures (A), (B), (C), (D), (E) respectively. 
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Digital visualization of hedgehogs 




Figure 5.14 : Neighboring dynamics. Figure (A) shows the rotatory motion 
of iterates in the Siegel compactum, while outer orbits are shown in (B) and (C). 

6 Conclusions 

In author's opinion, well-drawn figures of local invariant sets necessarily want 
customized approaches: equi-potential approach is looking like the best per- 
forming today and is could be even optimized when topologically equivalent 
model are applied; thus a general method to be applied for displaying a 
given local invariant set, independent on the input map, seems a utopia. 
But, being open to the possibilities offered by Science, one optimistically 
says 'up to now /' 

Alcssandro Rosa 
zandor_zz@yahoo . it 
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